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In this paper an attempt is made to find an interface of the perturbative BFKL Pomeron 
with the non-perturbative Pomeron originating from non-perturbative QCD phenomena 
such as QCD instantons and/or scale anomaly. The main idea is that the non-perturbative 
Pomeron involves a large scale (Mo m 2 GeV), which is larger than the scale from which 
perturbative QCD is applicable. One key result is that even for processes involving a large 
hard scale ( such as DIS ) the low x behavior is determined by an effective Pomeron with 
an intercept having an essential non-perturbative QCD contribution. 

I. INTRODUCTION 

In this paper we present an attempt to solve one of the most challenging problems of QCD: the interface 
between long distance (non-perturbative, soft) interaction and the short distance one (perturbative, hard) 
at high energy. The short distance interaction at high energy is controlled by perturbative QCD and can 
be described by the BFKL equation [1]. We call the hard Pomeron a solution to this equation, in spite of 
the fact that this solution is not a Regge pole. The region of long distances we describe by introducing the 
phcnomcnological soft Pomeron - a Regge pole with the intercept close to 1 [2] . Recently, progress has been 
achieved in a theoretical understanding of the structure of the soft Pomeron [3-9] , which provides hope for 
developing a first approach to understand the interface of these two Pomerons: the soft Pomeron and the 
BFKL one. 

This approach is based on a new hierarchy of the scales in QCD suggested in Ref. [3] and depicted in 
Fig. 1. 
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FIG. 1. Scales of QCD accordingly to Ref. [3]. 
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The key idea underlying this picture is the large scale for the soft Pomeron (qs — Mq ~ 2 GeV) which 
is closely related to the violation of the scale anomaly of QCD. The value of Mo is so large that the QCD 
coupling constant as (Mo) <C 1, and it can be considered to be a small parameter. From Fig. 1 we see that 
such sequence of scales makes possible a theoretical approach to the problem of interface of soft and hard 
Pomerons. For all virtualities (transverse momenta) larger than the separation scale, q > qn ~ 1 GeV , 
where as(q 2 ) < 1, we can apply the perturbative QCD approach. 

We formulate the problem of the interface between soft and hard Pomeron as a problem of a joint descrip- 
tion of the perturbative BFKL Pomeron with the soft Pomeron based on the non-perturbative contributions 
with sufficiently high scale (see Fig. 1). Many of the non-perturbative contributions originate from very long 
distances of the order of the confinement scale (see Fig. 1) but we assume that all of them are irrelevant to 
the structure of the soft Pomeron [3] . 

We also assume that the energies are so high that the saturation scale [10,11] is much larger than the 
scale of the anomaly scale, (or the scale of the soft Pomeron) , qs = Mq w 2 GeV and the separation scale 
qn which characterizes that value of the virtuality (transverse momentum) from which we can start using 
perturbative QCD. Our notation is transparent since qs gives the scale of the soft Pomeron, while for q > qn 
we can use the BFKL Pomeron to describe the high energy scattering in pQCD. 

In the next section we describe the general formalism of our approach which leads to a new Regge pole as 
a result of an admixture of soft and hard Pomerons with the trajectory as-H{t) = 1 + + a 's-H 

Section 3 is devoted entirely to the properties of the Green's function of the BFKL Pomeron with the running 
QCD coupling constant. As it has been shown in many papers [10,12-18] the QCD running coupling constant 
changes the character of the singularity in the BFKL Pomeron which becomes an essential singularity at 
u> — > instead of a pole at u> — ujl- In this section we present a detailed analysis of the Green's function 
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for the BFKL equation which is similar to that of paper [18] which appeared after we had completed our 
calculations. 

In section 4 we calculate the trajectory of the soft Pomeron which leads to the phenomenological Pomeron 
(ctphenom.p = 1 + 0.08 + 0.25 GeV~ 2 t [2]) and which appears as a result of the admixture of the BFKL 
and soft Pomerons in our approach. The deep inelastic scattering in the region of rather low virtualities of 
photon is discussed in section 5 where the values of virtualities for the hard processes are specified using our 
approach for admixture of "soft" and "hard" Pomerons. 

II. GENERAL APPROACH 

In this section we outline the general approach for dealing with the interface of soft and hard Pomerons. 
To do this we need to specify our description of the soft Pomeron based on general ideas of Ref. [3] . We next 
review the key feature of BFKL Pomeron in the LLA, elaborate its Greens' function in the Airy model, and 
finally discuss unifying soft and hard Pomerons. 

A. Soft Pomeron 

As a working model for the soft Pomeron we shall adopt the QCD instanton approach [4]. This is in fact 
not necessary since we only need to use general properties of the soft Pomeron exchange depicted in Fig. 2. 




Our approach: 

K(q 2 ,q' 2 ) = 

A S <l>(q 2 ) <l>(q' 2 ) ; 

^(q 2 )=e 73 ; 

q 2 = M 2 = 4GeV 2 ; 



FIG. 2. The main properties of the soft Pomeron originating from the QCD instanton approach [4] 

As shown in Ref. [4], the structure of the soft Pomeron in the instanton approach formally reduces to 
summing the ladder diagrams of Fig. 2. These ladder diagrams describe the i-channel gluon that propagates 
through the space inducing a chain of instanton transitions between different vacua and producing a num- 
ber of gluons in each transition. The sum of these diagrams leads to power-like (Regge-type) asymptotic 
amplitude, namely the scattering amplitude A(s,t) is asymptotically equal to 

A(s,t) = g(t)s a ^(t) ; (1) 
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+ a' s t. It is well known (see Ref. [19]) that both A5 and a' s can be found as 



with a S p(t) = 1 + A s + a' s t 
solutions of the Bethe-Salpeter type equation resulting from these ladder graphs. In this paper, we shall 
concentrate on the forward scattering limit, t = 0, where the kernel K(q 2 ,q' 2 ) is a function of transverse 
momenta, q and q', along adjacent rungs of the ladder in Fig. 2. We assume that this kernel has the general 
factorized properties, namely K(q 2 , q' 2 ) — As (j)(q 2 ) 4>(q' 2 ) which follow from the instanton approach of Ref. 
[4] but this could be more general. 



B. BFKL Pomeron 

The BFKL Pomeron gives the high energy asymptotic behavior for the scattering amplitude in the so- 
called leading log(l/x) approximation (LLA) of perturbative QCD [1]. In the LLA, for each set of Feynman 
diagrams of the order of a™, one keeps only leading terms proportional to (a s ln(l/z))™, and sum over these 
terms. In LLA this sum formally reduces to a summation of "gluon ladder" diagrams, (see Fig. 3). 
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FIG. 3. The BFKL Pomeron in LLA. 



However, due to the specific property of the vertex [1], the sum of these "ladder" diagrams leads to an 
amplitude A(s,t) oc s Ah / ^/ln(s), with Ah oc a s . The presence of the additional l/^/ln(s) factor reflects 
the fact that the BFKL Pomeron is actually not a Regge pole. 

For further discussion, it is more convenient to work in the complex angular momentum plane, I = 1 + lj, 
and to restrict ourselves to the forward limit t = 0. In particular, we need to find the Green's function 1 
G w (qf,qi) for the BFKL equation which satisfies the following equation (see Refs. [15,16,18] for details and 
Fig. 4 for graphical form of the equation) 



uG u (q f ,qi) = 5 {2) (q f -q t ) + J d 2 q K{q f ,q)G u {q ',$) . 



(2) 



x We shall also denote this Green's function as to emphasize that it is responsible for the "hard" perturbative 
QCD contribution. 



4 





FIG. 4. The BFKL equation for the Green function 
The kernel K(qj, <f ) in Eq. (2) can be written, after integration over the azimuthal angle, in the form 



K(q f ,<f) = a,(qj) K 



It can be shown that functions 

satisfy a generalized eigenvalue equation 



(3) 



(4) 



(5) 



with u)(f) as eigenvalues. Here, r = ln(q 2 /A 2 ) — ^\na s (q 2 ), m = ln(gfj/A 2 ), and ui(f) in leading order 
of pQCD is equal to 



«(/) = as(rH)Nc {2^(1) - - f) - + f)} 

7T Z Z 



In this paper we will mostly use the expansion of u>(f) at small values of /, namely, 

cj(f) = uj l (1 + Df + 0(/ 4 )) . 



(6) 



(7) 



It should also be stressed that in Eq. (6) and Eq. (7) the normalization point for the QCD running coupling 
constant a s was chosen as the separation scale, namely a s (q 2 H ). 

The solution to the inhomogeneous Eq. (2) can be found by expanding G u (rf, Ti) in terms of the complete 
set of eigenf unctions, Eq. (4), 



a + ioo 



df 



G u {q,qi) = G u (r,n) = / — x(w, /,/•*) <Pf{q ) = / — A x(u,f,n) 



a + ioo 



df 



2ni 



G(Y;r f ,n) 



duj 
Z~k\ 



(8) 
(9) 



where Y = ln(s/g/ 



5 



Using Eq. (8) one obtains the following equation for x( w >/> r i) [16] 

— = r H uj(f) x(u,f,n) + ne 1 . (10) 

Eq. (10) is a first order differential equation and can be solved analytically. The solution to Eq. (10) can be 
written as the sum of two terms, [16]: 

Xi(oJ,f,ri) = — df e J f , (11) 

2UJ J; 

and 

X2(u,f,n) =~2^ d f e 1 > ( 12 ) 

where the integration path in Eq. (11) is a ray Ci, from f to infinity, lying between the polar angles 7r/3 and 
7r/2 and the corresponding path for Eq. (12) is another ray Ci from f to minus infinity, lying between the 
polar angles — 7r/3 and —ir/2. For G w (r/,rj), we obtain 

rf+ioo if 

G u {r f ,n) = -± e* r f (xi(u,f,n) + X2(u,f,r i )) . (13) 

Jf-ioo 27Tt 

It is worthwhile mentioning that the difference between Eq. (11) and Eq. (12) is a solution to the homogeneous 
BFKL equation with the running QCD coupling constant, (Eq. (10) without the inhomogeneous term). The 
fact that x(w, /, r,) = xi(w,/, r,) + X2(^, /, r i) corresponds to having chosen the integration constant to 
Eq. (10) in order to ensure that G w {rf,Ti) is real. 



C. The BFKL Green function in the Airy model 

To clarify the main properties of the Green's function we consider the expansion of Eq. (7) which is justified 
at small values of / and /'. We call this approach the Airy model, and it was first done in Ref. [18], our 
functions will be very close to the Airy functions. Using new variables /+ = /' + / and /~ = /' — /, we 
can reduce Eq. (13) to the form 

G^(r f ,n) =g{ U ,r f ,r t ) = -p- ( f df~ e^J-^u) + f df - e ^J-,r f ,rA (w) 

4^ Jf_ ioo 2iri \J Cl J C2 J 

*( W ;/+,/-;r / ,r i ) = -^/- + + {1 + f [\i.f + f + \ (.rf}}; (15) 

Integration paths C\ and C2 in Eq. (14) are defined in the same way as in previous section, with the initial 
point shifted to /~ = 0. 

Since the /+-integral is Gaussian, we have 



IQiruiLUL th D \J Cl 



Co 



g{u,r f , n ) = W- ! + / -= e > « ^ ' »« f~' " 



(16) 

1/3 



with the Co — 1"+^ an d Sr = rf — ri . We can further change the variable to: v = I du>lt h ) f 

2 ^ ' 

and obtain for our Green's function: 
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M&Q; (17) 

(18) 
(19) 
(20) 

Note, that the right dimension of solution of Eq. (2) must be g(u>,rf,ri) / (q/Qi) , but this fact is not 
important for our later consideration, so we will proceed with the dimensionless function Eq. (17). Once 
g{u,r f is known, the asymptotic behavior of our Green's function G(Y;rf,ri) can be obtained from 
Eq. (9). However, in dealing with the general expression ( see Eq. (17) - Eq. (20) ) we have to take into 
account that the limit depends on the virtualities rj and ri and, in particular, the value Co. Further, 
ultimately, we need also to take into account the soft Pomeron intercept. For different relations between 
these two values, different region in the uj integral can dominate, leading to different asymptotic behavior 
for the Green's function. We start the search for the resulting singularities by considering various situations 
where different regions of ui dominate the integral in Eq. (9). 

We first consider the case when uj » Co. In this region both £ ( Eq. (19)) and ( (Eq. (20)) are large, and 
therefore the main contribution in this integral comes from the region of small v. Neglecting the v z term in 
the exponent of the integrand, we obtain 

A(€;C) oc r^e-vi-i? (21) 
Jo V v 

We will see later, that omitted e v ^ 3 term is indeed small. For such an integral the integration over two 
different contours in Eq. (13) give the same answer, since the integrand has no singularities in the right 
semi plane and both integrals could be taken along the real axis. The resulting integral can be evaluated 
explicitly, leading to: 

Mn,r f ) cx 2if- 1/2 ^(fr/2)^/^j (^) 1/2 (22) 

where K u (z) is the modified (McDonald) Bessel function. Since Eq. (21) does not depend on the sign of Sr, 
we can therefore write the answer as 



g(co,r fl r t ) = Gi(r fl r t ) = *J — - (e- 2 ^^£ 6(5 r) + e 2 ^^^ 9(-S r)) . (23) 

Instead of integrating over Sr, we could simply integrate over the following function for positive values of Sr: 



Gi(r f ,n) = « J e -2(¥)^ , (24) 

r 1 ' 2 V lulolth D 



where f = Ti + 2 r ' s — ^ L J H ■ In this region, Sr and v are small, Sr <~ yj J 3 ^- << 1 and v cx jfj? ~ 

1/3 

( 3§7u) <<; 1- Therefore, we can trust this solution so long that the following condition is satisfied: 

- < 1 • (25) 
to 

Indeed, even when ^ = 1, the v is still small, v <~ (-gj) 1 ^ 3 << 1 • The magnitude of the omitted term in 

j p a _ 

96 rf + r-j 

our solution is also small, namely, e " ( 2 ' , and we have for the Green's function in this region: 
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This equation gives the asymptotic answer for the Green's function for the case u> » Co, in the so called 
diffusion approximation for the BFKL kernel. The Mcllin transform of Eq. (26) gives: 



G\Y-r f ,r i ) = x l^- v ^e-^ 



(27) 



( see calculations in Appendix A ). 

There is another contribution in the limit of u » u) coming from the large v region. Indeed, it is easy 
to see from Eq. (18), that there is another saddle point, 



y SP 



£ * 0. 



In the region £ >> 1 we have 



vsp = V£ = r) — ] W {ulThD/A) s >> 1. 
This saddle point corresponds to the second crossing point of Fig. 5. This is a solution for the u > Co. 



(28) 



(29) 



1/ocP co » co 




Y/r 1 / 2 



CO - C0« CO 

1 2 3 4 5 6 7 CO/ CO 

FIG. 5. The positions of the saddle points in lj for G(y;rf,ff) for lo > u>. 

However, it turns out that in this regime, fg p » 1 while fg p is small (fg P <C 1). It follows that we can 
not use the diffusion approach for the BFKL equation to obtain an estimate of this contribution. In this 
case Eq. (17)-Eq. (18) are no longer valid, and we have to take the full kernel of Eq. (6) into consideration. 
The integration over this full kernel is performed along the real axis and instead of two solutions Eq. (11) 
and Eq. (12) now we have one. Going back to the variables f + and /~ and considering /+ / /~ << 1 , one 
can obtain the following equation for the saddle point: 



r f ~ 



2r H a s (r H )N c 1 



1 - /■ 



= 



(30) 



SP 



Eq. (30) stems from the fact that in the region of large w, f is rather large. For large / we can replace 



co(f-) byw(f-) 



a a (r H )N c 



-. Thus we have 



(31) 



We obtain for the Green's function in this region of uo: 



Gl(n,r f ) ~ Veoe-^- 2 ^^" Kg) (32) 



For the w-saddle point, we have 



which leads to 



TTLOgp LO$P 



,'2a s (r H )N c r H 2a s (r H )N c r H \ 
"SP = ( ^ y-r- 2 ) ■ (34) 



Using Eq. (34), we arrive at 



/ 2 Qs (r g )JV e r g 2 as (r H )JV c r B 

G 2 ^;^) ~ e 2 V - ^ _ (35) 



This contribution has been missed in Refs. [16,18]. Nevertheless, comparing Green's functions of Eq. (26) 
and Eq. (32), we see, in the limit of small Co, the contribution of Eq. (26) is larger than the contribution of 
Eq. (32): 

^4 ~ e"^ . (36) 

We next consider the limit where u) is close to Co: uo = Co + 5co, i.e., 5lo / Co « 1 . In this case £ is small, 
£ << 1, indeed, £ <~ oc 4£ << 1 , and for small £ we have v ~ 1. We also take 5r in the limit 

<5r << (rjj D) 1 / 3 and, therefore, we neglect the £ 2 jv term in our equation. We now need to consider both 
solutions of Eq. (10) together, sum of Eq. (11) and Eq. (12). Due to the integration along the rays tt/3 and 
— 7r/3 we have : 

- £ v e 



and finally we obtain ( see calculations in Appendix A ): 



e-^ e (37) 



m ( 4a; (^) 2 „ , 2 r(|) 



The Mcllin transform of Eq. (38), see Appendix A, gives: 



The region of applicability of this saddle point analysis stems from condition that /+ and /~ should be 
small in order to use the expansion of Eq. (7) and from lo sp < Co. It is easy to see that these conditions 
restrict the value of rapidity Y [13-16,18] 

r 2 

< Y < — f — . (40) 
lol r H 

Now, if we consider together the asymptotic behavior of the Green's function given by Eq. (27) and 
Eq. (39), we see that in a large range of rapidity, given by Eq. (40), the contribution of Eq. (39) is much 
larger than contribution Eq. (27). Therefore, the main contribution to the Green's function gives the region 
of such lo, where uo <~ Co. 

In the region of small lo < Co where £ is negative (see Eq. (19) ), the saddle point of v is equal to 
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v S p = 




Performing calculations for the saddle point we arrive at 

j2uj L r H 

usp = ± (1 -*) y D3Y • (42) 

The Green's function is [18] 

G(Y;r f ,r f ) ~ e 2 cos . (43) 

It turns out that restrictions fg P < 1 and fg P < 1 can be satisfied, but only within numerical accuracy, 
namely, fg P < 1/D. We will see later that D is rather large. 



D. Soft to BFKL Pomeron transition 



The main idea of the paper is shown in Fig. 6. Generally speaking the sum of diagrams in Fig. 6 can be 
found as a solution to the following equation: 

Gl- H (r f ,n) = D{w;r f ,n) + D (g) (K s + K H ) (g) G^" . (44) 

where sign (^) stands for all needed integrations and D(u>;rf,Ti) denotes the Green's function of two gluon 
exchange. Using the factorization properties of soft Pomeron kernel we can solve this equation, the derivation 
of the solution is given in Appendix B, and it reads : 

G S ^ H (r f ,r t ) = G*(r f ,n) + - _ — j — j X-'/(r')W(^^") ' ^ 

where 

G^(r f ,n) = As J dr"<j>(r") G^{r f ,r") J rfr'^(r') G^(r',r 4 ) . (46) 

It is easy to see from Eq. (46), that when Ks = , Eq. (45) describes only the "hard" Pomeron. In the 
case, where there is no pQCD emission, when G H = D , Eq. (45) has a form: 



/ dr" 0(r") D(Lj;r f ,n) J dr' cf>(r') D(u;r f ,n 



°s < r " r <> - D ^ ri) + As i- Aw-^y^o^;^^;)^) • (47) 

Since D(u;rf,ri) — 5(rf — ri)/u) , Eq. (47) leads to 

r<s-H (r r s S ( r f - r i) , 1 K s(r f ,n) 

G u {rf,n) = + -— -r— , (48) 

or, in other words, we only have the soft Pomeron in this case. 



10 



FIG. 6. The diagrams and the graphic form of the equation for soft to BFKL Pomeron transition. 
Eq. (45) has poles which are the zeros of the denominator, namely, 

1 - A s J dr' J dr" (f>(r') G^{r' >") <j){r") = 0. (49) 

For uj » u> the solution of Eq. (49) is (see Appendix C ) 



96rUu s + A s ) 



lu = A s _h = A s + ci s + ng 2 t \ , . (50) 



Therefore, the intercept of the resulting pole is simply As-h = As for uig — * . This result is expected 
since for ujs = 0, there is no perturbative emission of gluons. It should be stressed, that this answer is 
obtained with the restriction, 
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w > lo s , (51) 

see also Eq. (25). 

It is easy to see, that this restriction is satisfied when u>s is small. This happens in two cases: 

• when as is small. In this case w is also small and, therefore, from Eq. (51) we have: 

A s >> S ~ w • (52) 

We can make a simple numerical estimate: for as ~ 0.1 we have lol ~ 0.18 (d>s w 0.1) (see next 
section), so taking w = As-h ~ 0.2 we obtain As ~ 0.1, where As > ws ~ w. So, one can see that 
restriction Eq. (52) is satisfied for as < 0.1; 

• when Cos is small due to possible large rs, namely, 

rs > -r • (53) 

A S -H 

We next consider another region of to. We assume that the position of the resulting pole to = A$-h may 
be close to ujs , As-h = los + 5lo , where for Slo there are two asymptotic limits. The first is the small Slo , 
%^ « 1 , and the second for the 5lo of the order of Cos , %^ ~ 1 • The solution of Eq. (49) in this case is 

»- A »-- fi « + *'-^(^)'- < 54 > 

Fig. 7 shows the numerical solution to Eq. (54) with rs — r™ at (see later Eq. (58)), and with As = 0.1 . 
We obtain that for A, = 0.1, ~ 0.3 - 0.8. 




12 



FIG. 7. Solutions for As-h and d>s as functions of as at A s = 0.1. 



We can trust this solution only if to is larger than ujs , and this condition can be translated in the following 
inequality 

2A%r s ' ! 



**< af^laWj ' (55) 

or in the restriction on the value of ws 

_i/3 9V3 rfll 

«5<A fl ^^-^. (56) 

Remembering that Eq. (54) is valid only if ws > As we see that Eq. (55) is a solution for rs which satisfies 
the inequality 

1/3 2 V3 lYl) 

A S < us < As-^^AlZ. (57) 

The solution of Eq. (57) defines the region of rs for which we can trust the approximate equation Eq. (54). 

In Fig. 8, we present the solution for rs to Eq. (57) as a function of a$ at fixed As = 0.1. In Fig. 9, 
the same solution is plotted versus ^ for D = 3 . 




123456789 10 

co L /A s 



FIG. 9. Permitted rs as function of at D = 3 . 



The matching between Eq. (54) and Eq. (50) occurs at rs = r™ at , for which we have 



Du>s A s 2A|r s / T(i 

96r|(% + As) u 2 s D 2 \3 2 / 3 



, S -h = A s + u s + 7T7 - T7T ^ rT = *f [ ^ ] • (58) 



A numerical solution of this equation is presented in Fig. 10. This solution is the function of as with 
As = 0.1 . To justify this procedure, we have to show that, for rs = r™ at both solutions, Eq. (54) and 
Eq. (50), are valid. It is easy to see from Fig. 8 that r™ a * calculated from Eq. (58) for any value of as > 
is in region where Eq. (54) has solution. From Fig. 7 we also obtain, that for given rs the ratio 4^- is 
4j- <~ 0.3 — 0.8 , and these values satisfy the condition of Eq. (A14), which defines the applicability of this 
solution. The solution Eq. (54) is valid for rs = r r g lat . We next consider the solution given by Eq. (50). 
The region of applicability of this solution is defined by Eq. (53): 

^ UJLrH 

rs > TT-r (59) 

^ As-ff 

The boundary value of rs, which is rs = 2 ^ s rff H , is also presented in the Fig. 10. Considering this picture, 
we see, that this line lies under the r™ at defined by Eq. (58). Remembering, that the ratio 4^ ~ 0.3 — 0.8 
satisfies the Eq. (25), we obtain, that solution Eq. (50) is also valid for the given values of r™ at . 
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0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 

a s 

Matching value r™ at and boundary value r$ as function of as at As = 0.1. 

FIG. 10. 

So, finally, we have two solutions, the first is Eq. (50) which is useful for r$ that lies between the lines of 
Fig. 10, and second is Eq. (54), which we use when rs lies between the line for r s 1,at from Fig. 10 and upper 
bound of Fig. 8. In the following we will denote all these solutions by the As-h- 



III. NLO BFKL 

As have been discussed the position of the resulting singularity (pole) crucially depends on the value of Co 
or uj l . In this section we discuss the values of ojl and D in the NLO BFKL equation [20,21]: 

lu = uj l (l + Df). (60) 

NLO BFKL equation leads to the following values of wi and D (see Refs. [16,21,22]): 

uj L ^a s (gg) (2.772 - 18.3 a s ), 

and 

D = a s (q%) (16.828 - 322 a s ) , (61) 

These values have a serious pathology, as for as > 0.152,(g§ < 30GeV 2 ), ljl becomes negative and 
for as ~ 0.05 , (q 2 < 3 x 10 6 GeV 2 ) D also changes sign. A solution to this problem was proposed by 
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G.P.Salam [23] and M.Ciafaloni, D.Colferai and G.P.Salam [17]. They suggested a scheme to improve 
(resummation) the NLO BFKL kernel. Their resummation is based on the idea that the resumed NLO 
BFKL kernel should reproduce the NLO DGLAP kernel in the region of large photon virtualitics. Indeed, 
when we start to consider the NLO terms 

a s [ a s In — ] (62) 

and change the energy scale in the DGLAP evolution equation, namely, s = k\ {k\ » fc| ) to the energy 
scale s = k\ k 2 typical for the BFKL approach, the double logarithmic terms arise: 



ln2 d> K ln 7i ln 7^-) • (63) 



£2\ / ^2 

^2 / V * 11 ^2 ^2 

These terms, formally subleading, nevertheless may give significant corrections to the NLO kernel. It was 
shown [23,17] that these terms are responsible for the difficulties with the NLO kernel. Therefore, the 
proposed resummation of these double logarithmic terms was done so, that the resummed kernel reproduces 
a correct DGLAP limit for lu ± 1/2 and for a scale sq = kf . Using this method of resummation we can 
avoid the difficulties of the NLO BFKL kernel. This improved kernel, reproduces correctly the exact NLO 
BFKL kernel in the region of its applicability, but it turns out to be free from inconsistency of "pure" NLO 
BFKL kernel. From the proposed different schemes for the improved kernel, see [23] , we use scheme number 
four, as this scheme and also scheme number three has a reasonable behavior of second derivative of the 
kernel as a function of d s . As the additional resummation of the terms which are important for the second 
derivative of the kernel was made. 



So we have: 



Xo (7) - ../o - 1 /o , ../o + 1 /o -./o ../o , , 5 + 



1/2-7/2 1/2 + 7/2 l/2- 7 /2 + w/2 + a s B 1/2 + 7/2 + w/2 + a s B 



{ A (l/2-7/2 + 1/2 + 7/2 ~ 1/2 - 7/2 + ^/2 + a, B ~ 1/2 + 7 /2 + u/2 + a s b) ) - (64) 

Here 

X (7, w, a s ) = a s xo (7) + °? s xi (7) 
is usual BFKL kernel in the next-to-leading order, and 

B ~ 8 12iV c + 6iV3' (65) 

Solving the NLO BFKL equation with improved kernel 

w = X (7, a s ) , (67) 

we can calculate our NLO and D. The results are presented in Fig. 11 and Fig. 12, where we also plotted 
LO lul and D . These figures show that the resummed NLO BFKL kernel gives smaller values both for lul 
and for D. From Fig. 11 we obtain iD w 0.2 , which is too high to describe the experimental data for soft 
processes. On the other hand, namely, this value of the Pomeron intercept describes the inclusive data [29]. 
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FIG. 11. 
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FIG. 12. 



IV. HARD AND SOFT POMERONS INTERFACE IN DIS 



In this section we consider the deep inelastic processes in which we can observe our "effective" Pomeron 
considering the Q 2 - dependence of the cross section ( Q 2 is the photon virtuality) . The diagrams that 
contribute to DIS are shown in Fig. 13. 
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target 



target 



FIG. 13. Diagrams and graphic form of contributions to deep inelastic processes. 
The DIS total cross section can be written in the form 

<Jtoth*p) = [ dr f dnJ~<'(Q 2 ;rf) G DIS (Y,r f , n ) J N (r H , n ) , 



(68) 



where J 7 and J N are photon and nucleon impact factors, which can be easily calculated in the LO of 
perturbative QCD (see Refs. [24] and references therein). For example, 



d 2 r ± / dz\*{Q 2 ;r ± ,z)\ 2 (l - e^x ) 



■/ 7 *WV/) 



We introduce the Mellin transform of G D1 (Y,rf,ri) with respect to Y 



G ulh {Y,r f ,n) 



1 



a-\-i oo 



did 



If It Ja-i oo 2 7T i \ qf qi 
G® IS is sum of diagrams of Fig. 13 and, therefore, it is equal to G^~ H given by Eq. (45): 



G, 



DIS 



(TV rO - G H (r f rl + ± J dr "^ G% (r f ,r") f dr^(r>) G»{r\n) 
" K h %) + S 1 - A s J dr' J dr" 4>{r') G"{r', r") <j>(r") 



(69) 



(70) 



(71) 



The Green's function of Eq. (71) has a quite different asymptotic behavior than g (w, r/, r^) due to the 
contribution of the denominator. The zero of denominator generates a Regge pole ( resulting Pomeron ) in 
lu . Therefore, the asymptotic behavior of Green's function of the Eq. (71) depends on what is dominant in 
given kinematic region, either Gfj or this pole. 

First, we consider the region of rapidity, given by Eq. (40). Now, we have two conditions on As-h, which 
should be satisfied, so that the soft-hard pole contribution in Eq. (71) will dominate. First is the condition 
of the existence of soft pole solution in denominator of Eq. (71). This condition is given by Eq. (52)- 
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Eq. (57). Second , we need to check that this soft-hard contribution to the Green's function will give a larger 
contribution than either the contribution of Gfj (r/,r,) or of the numerator of the second term of Eq. (71). 
From the derivation of the Eq. (A15)-Eq. (39) is clear, that the asymptotic behavior of both functions in 
the leading order is given by e uY . So, all these conditions together define the following restriction for the 
kinematic region for resulting Pomeron : 



> Cj. (72) 



We have obtained this restriction before, see Eq. (52)-Eq. (53), and can say again, that this restriction is 
satisfied for the case of very small as , when as < 0.1 and therefore As > ~ u> , or for the case of 
large r ' ^ r/ , when Tf \ r ' 1 > £ L ra H ■ So, in both these cases asymptotic behavior of Green's function will be 
defined by resulting pole. Otherwise, we have the answer of Eq. (39): 

G < y; " • r " > = \l sJJir H n " ~ (73 > 

in the region of rapidity: 

r 2 

< Y < — I— . (74) 

Another possibility for the resulting Pomeron dominance is to have large values of rapidity. In this case 
the Eq. (73) is not valid anymore, and we have to use equations Eq. (41)-Eq. (43). So, using Eq. 42, we can 
claim that for: 

A S -h > losp (75) 

we obtain the constraint on the value of rapidity for which the resulting Pomeron governs the high energy 
asymptotic: 



Let's assume now, that we have resulting Pomeron regime, i.e. consider only the second term of Eq. (71). 
We want to rewrite our Green's function and estimate the resulting Pomeron contribution. For that, first of 
all, we have to expand the denominator of resulting Pomeron, Eq. 71, with respect to ( w — A S -h ) : 



1 - A s J dr' J dr"<P{r') G%{r',r") 0(r") w 1 - ^ 



Here we take 5u) and As-h from Eq. (50). Using this expansion, we rewrite our Green's function: 
MS,,, ... ^ _ ^ [ a+t0 ° d " ( s V 9{^r f ,n) 

where 



G ults {Y ,r f ,n) = / -fn __ ^ -21 (77) 

QfQi Ja-too 2ttz \q f q l J ui - A S -h 



g(w, r f ,n) = J dr"cb(r") G»(r f ,r") J dr' <j>(r') G»(r',n) 
Taking the contour integration we obtain the answer: 



As / s ^ As - H 



G DIS (Y,r f ,n) = ( — ) g(w = As-H,r f ,r i ) , (78) 
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or using the Green's function of Eq. (C3) we have: 

9/ Qi (oj l r H D ) 4/3 V 9/ * / 

So, from this consideration, we see, that we have the resulting, phenomenological, Pomeron contribution, 
determined by Eq. (79), in the following cases: 

• first, in the range of rapidity 

< Y < Tf 



for the small values of as : as < 0.1 ; 
• second, in the range of rapidity 



r 2 



for the large values of — 



Ti + rf 



< Y < 



rf + n lo l r H 

9 > A r5 

2 A S -H 



• third, for large rapidities: 



> 



DA% 



One can see from Eq. (77) that the soft-hard transition generates a Pomeron satisfying the initial condition 
for the moments of the deep inelastic structure function. Therefore, in the kinematic region, where the As-h 
contributes , the DGLAP evolution equations only affect the Q 2 dependence of the deep inelastic structure 
function, while the energy dependence is entirely determined by the soft Pomeron. Such a situation has been 
studied in details in Ref. [26] . In other kinematic regions we have a normal DGLAP evolution equation with 
the initial conditions which do not contain the soft-hard Pomeron contributions. 

It should also be stressed that at fixed virtuality Q 2 at high energy, asymptotic behavior is determined 
by the exchange of the soft Pomeron independent of the value of Q 2 . However, the coupling of the soft-hard 
Pomeron to the photon depends on the value of the photon virtuality. 

Estimating our function, Eq. (79), for the case of large Q 2 , we see that the contribution of the rightmost 
singularity ( soft-hard Pomeron) is small at large qj . This fact makes our life more complicated and we 
have to remember the infinite number of the poles that we have in Eq. (45) . Indeed, these poles accumulate 
to zero, as will be seen latter, and, therefore, we face the problem of an essential singularity at zero, which 
possibly may give a large contribution to the Green's function. We have to estimate the contribution of 
these secondary poles in G DIS ( Y ; , 77 ). This estimate is considered in Appendix D, and the answer is 

r niS (v . r , r«2-_V? y r iy> ((^& Ok=r s (M^ 0)r f =rsh=^ |(A)il^ ln _^_ 

° 1 ,rf,Ti) - AsV*D 5/6 V ( 7rn ) 4/3 

(80) 

The estimate of the value of this sum we may perform using the method of steepest descent. This estimate 

gives the eV ln ^ behavior for the sum, and this contribution is smaller, than the contribution from the right- 
most singularity, As-h- Therefore, in our calculation, we can only need consider the resulting Regge pole 
contribution and neglect the contribution of the secondary poles. 
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V. CONCLUSION 



This paper is a realization of old ideas to interface non-perturbative Pomeron and perturbative BFKL 
Pomeron [28] . The new ingredient in our approach is that the typical scale for the "soft" Pomeron is consid- 
ered [3,4] to be rather large ( M ss 2 GeV or rs — ln(Mg /A 2 ) rs 4.6 ) , larger than the values of transverse 
momenta (pt for which we can use perturbative BFKL approach (p t > lGeV,rH Rs ln(lGeV 2 /A 2 ) = 3.2). 
In other words , we assume that rs > rjy. 

We hoped that taking into account the running QCD coupling constant we would obtain the resulting 
singularity which would govern "soft" processes close to unity, as it appears in the experimental observations 
[2] . This hope is based on the fact that the BFKL Pomeron with running QCD coupling leads to a singularity 
at a angular momentum which is equal to unity. However, it turns out that even for running a s the resulting 
singularity is located to the right of uis = lj bfkl ^t-. Using the next-to-leading order BFKL estimates for 
uj bfkl rs 0.2 [20,21,23], we obtain ivs = 0.139 which is larger than the Pomeron intercept of the D-L 
phenomcnological Pomeron [2]. The estimate of the Pomeron intercept from the single inclusive production 
suggests a higher value of the Pomeron intercept [29] which is in better agreement with our prediction. 

The main result of this paper is concentrated in Eq. (50) and Eq. (54) which give the position of the 
resulting singularity for different ranges of a s and non-perturbative input As . We also discussed the interplay 
between the BFKL Pomeron and the resulting Pomeron in deep inelastic processes which provide a possibility 
to measure the hard BFKL Pomeron at large virtualities of photon, and to investigate the interface between 
the "soft" and the "hard" Pomerons for small values of the photon virtualities. To our surprise our result 
shows (see Eq. (79) ) that even at large photon virtualities the contribution of the resulting singularity 
dominantes at high energy (low x ). In other words we cannot consider DIS at low i as a typical hard 
process where we can safely use perturbative QCD. The admixture of non-perturbative interaction is so 
essential that it changes the effective Pomeron intercept. 

It is interesting to compare our model with the model introduced in Ref. [18], where two Pomeron model 
was also introduced. We have no ideological disagreements with this model but we differ from it in the way 
a transition from soft to hard Pomerons occurrs and in the origin of soft Pomeron. We assumed that the 
existence of the soft Pomeron is related to typical non-perturbative phenomena e.g. QCD instantons and/or 
the scale anomaly in QCD, while in Ref. [18] M. Ciafaloni, M. Taiuti and A. H. Mueller study a possibility 
that the non-perturbative Pomeron would be generated by the behavior of the QCD coupling at low values 
of transverse momenta. Our "soft" Pomeron has normal Gribov diffusion in impact parameter space which 
we admixture with diffusion in transverse momenta typical for the BFKL Pomeron. The freezing of the 
QCD coupling leads to a new picture which was suggested in Ref. [27] and which is certainly an alternative 
approach to ours. 

We would like to stress that in this paper we consider only the interface of two Pomerons: the soft Pomeron 
and the BFKL Pomeron. The reality could be much richer and could possibly lead to another scenario of 
interface the long distance interaction with the short distance QCD physics. As was pointed out in Ref. 
[30], the soft Pomeron could be suppressed if it is built from QCD instantons, [4], in the high parton density 
environment. In other words, if the interface between long and short distances interaction passes the stage 
of high parton density QCD, the soft Pomeron does not contribute to the high energy asymptotic behavior 
of the scattering amplitude. 

The main result of this paper we consider to be the observation that the non-perturbative correction could 
be essential for so called "hard" process at low x and large scale of hardness ( like DIS at low x). 
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APPENDIX A 



We consider the mellin transform over u> in Eq. (26) with the contour of integration defined in a usual 
way: to the right of all singularities of the integrand. First consider the limit Y — > 0: 



2r, n~ r +lo ° duj e~ 5r ^ 



G\Q;,r f ,n) = — l — J— / — = . (Al) 

n+r f \ ljD J Q _ ioo 2m s/w - w 

Closing the contour of integration to the left over the cut from the — oo to Q with the variables change: 

lo - lo = e l7T 5lo (A2) 

on the upper edge of the cut, and 

lo - lo = e~ iir 6u (A3) 
on the lower edge of the cut. In so doing, the integral Eq. (Al) turns into a sum: 



or 



G\0;r f , ri ) = 4n r - ^ f ^ , (A5) 

^ } 1 n +r f J Q 2tt n +r f J x 2ir V ' 

which leads to (after changing variable Slo — > — Slu in the second term of this expression) 

V f ' ' n+r f 2tt n+rf y ' K > 

Coming back to Eq. (Al) with non zero Y, we have 



(y;r/ ' ri) = ^7 L ioo 2^1 7^ ■ (A7) 

This integral can be evaluated using the method of steepest decent, with the the saddle point in lo determined 
by: 

^-TTTFT^ ' (A8) 

In the limit lo » Gj , Eq. (A8) leads to the following saddle point value: 

5rV 1 1 



USF= It (A9) 



and to the asymptotic behavior of the Green's function: 



^;r/,rO = V^^e-^, (A10) 

which is valid for ujsp » lo, or lo ~ 1/Y » lo, or Y « A. 

Now we consider Eq. (37). Here there are complex conjugates functions and, therefore, the answer is the 
sum of the real parts of these functions: 

g(LO,r f ,r f ) cx 2 cos(tt/6) / -= e~ iv cos(7r/3) " t cos( £ z/ sin( J)) + 
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2sin(7r/6) r ^ e-^coe^/s)- 4 sin (^ sin(J)) . (All) 
Jo V ^ 

Because £ and <5w are small, while f is the order of 1, we can expand our function with respect to £ v- term. 
Integrating this expansion over v we obtain the following result for the first two orders of £ : 



Taking into account expression for £: 



we obtain finally: 



3 2 / 3 \Du L r H J 8w 2 v ; 3 4 / 3 



This answer is also valid for the case when Suo /Co <~ 1 and £ ~ 1. In this case, the series of Eq. (A12) 
will be convergent due to the in e~^ v expansion. Indeed, we see from Eq. (A12), taking £ = 1, that 

r(§) 1 

6 a — « 1 , 



T(i) 8 3 2 / 3 80 

is small. The following condition for £ arises in this case: £ < 3, and we obtain the simple estimate for 
the region of oj where the series is still convergent and the given solution is still valid : 

% < 3 , (A14) 

where we still have < 2 . Now, if we go back to the Eq. (26), and taking u = Co, we see, that the 
correcting term to Eq. (26) is still small, 



1 + V, <<1 

Therefore, in the case Slo /Co <~ 1 and ^~ 1 we obtained overlapping solutions for Green's function, Eq. (26) 
and Eq. (A12), which match each other. So, we now know the form of our Green's function in the whole 
region of positive £. 

We also find here the G(Y; r/ , rj ) function. Performing Mellin transform of Eq. (All) and integrating 
over Suj in the first approximation in t;v we have : 

G(Y;rf,n) cx / — -e V »«l r H J 2*2/3 w ' / _ e 3 ) (A15) 

J ice 2tti Jo V" 

or 

G(Y;rf,ri) cx e / — — -e V »»l r H ) 2^/3 / — = e 3 i (A16) 



or 



Gfyjr^rO cce^f^SiY-, (j^)^ ^ -(./a) ) . (A17) 
The final answer, obtained for this Green's function in this asymptotic limit, is: 



G (»/,r0= ^/;; rg ^ -^<^^. (A18) 
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APPENDIX B 

The Eq. (44) has the following form : 

G S ~ H (r }l n) = D(Lj;r f ,n)+ J J D(w;r f ,r a ) (K s \r a ,r b ) + K H \r a ,r b )) G S - H (r b ,n) d 2 r a d 2 r b . (Bl) 
We are searching the solution of Eq. (Bl) in the form 

G s - H (r f ,n)= J G H (r f ,r a )f(r a ,n)d 2 r a . (B2) 
Inserting Eq. (B2) in Eq. (Bl) we obtain: 

J G H (r f ,r a )f(r a ,n)dr a = ^ + ^ <f>(r f ) J J <j>(r a )G H (r a ,r b )f(r b ,n) d 2 r a d 2 r b + (B3) 

+ JJ jD(cu;r f ,r a )K H (r a ,r b )G S - H (r b ,r c )f(r c ,n)d 2 r a d 2 r b d 2 r c , 

where we used that K s (rf,ri) — As <f>(rf) <t>{r%) and D(u>; r/,r,) = S ^ rf a) ~ r '' > ■ 
And because (see Eq. (2)) 

G H (r f ,n) = D(Lj;r f ,n) + J J D(u;r f ,r a ) K H (r a ,r b ) G H \r b , n ) d 2 r a d 2 r b (B4) 

we have: 

J G H {r f ,r a )f{r a , n )dr a = ^ 'J n) + ^j- (f>(r f ) J J <P{r a )G H {r ai r b )f{r b ,r t )d 2 r a d 2 r b + (B5) 

+ J (G H {r f ,r a ) - ra) ) f(r a ,n)d 2 r a . 

This equation leads to 

f(r f ,n) = S 2 (r f - n) + A s 0(r/) J J 4>(r a )G H (r a ,r b )f(r b ,n) d 2 r a d 2 r b . (B6) 

Multiplying both parts of Eq. (B5) on the integral J G H (rf,r a )<p(r a ) d 2 r a and integrating over rf we 
obtain: 

J J <P{r a )G H {r a ,r b ).f{r b ,r i )d 2 r a d 2 r b = j G H (r a ,n)^(r a ) d 2 r a + (B7) 



(B8) 



+ A S J J <j>(r a )G H (r a ,r b )(f>(r b )d 2 r a d 2 r b J J <j>(r a ) G H (r a ,r b ) f(r b , n) d 2 r a d 2 r b . 
Eq. (B7) gives: 

( \nH< \fi \fi .ft J G H {r a ,r t )(j)(r a )d 2 r a 

<{r a )G (r a ,r b )f{r b ,ri)d r a d r b = — - — . 

I - A s J J <j>(r a ) G (r a , r b ) 4>(r b ) a 2 r a a 2 r b 

Now, putting back Eq. (B8) into Eq. (B6) we obtain: 

f(r f r) S 2 (r f - r) + ^ <f>(r f ) J G" (r a ,n)cf>(r a ) d 2 r a 

m, t) [ f l) u 1 - A S J J 4>(ra)G*{r ai r b )^r b )d 2 r a d 2 r b ' ^ 

And, finally, we obtain for the solution of Eq. (Bl) 

r S-H (r r N r H (r \ | a I G H (r h r a )<j>{r a )d?r a J G H {r^n)^) d 2 r a 

G (r f ,n) = G (r f ,n) + A s j - - ; ; _ _ ^ - ^ - . (BIO) 
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APPENDIX C 



We search the poles which are the zeroes of denominator 

1 - A s J dr' J dr"cp{r')G^{r',r")<P{r") = 0. (CI) 

To calculate Gfj (r', r") in this equation, we use the result of Eq. (26) and Eq. (38) for different regions of 
co. For oj » lj we can use the following asymptotic expression: 



r 1 / 2 V u)LU L r H D 

while for u> ^ a), £ << 1, « 1 , as well as for £ <~ 1 , <~ 1 we can use the asymptotic answer 
for Green's function given by Eq. (38): 

g 2 (r> r") - r - 4 / 1 r 4 " v /6 f Eiii r 2/3 ( r ^) 2 ^ nil (C3 ) 

GJr ' r) - 4 i uu L r H D \u L r H D) \Du L r H ) 8^ ^ 3^ J' ^ 
Let us now consider Eq. (CI) using Eq. (C2) for (r' , r")\ 

1 - A s J dr' J dr" 4>{r') Gj, (/,/') 4>{r") = 0. (C4) 



Using new variables 



in Eq. (C4), we obtain: 



x = — t, — . y = — o — . ( C5 ) 



-^e^,fc^ Hay r tee-y^e-*"*™" . (C6) 
w - w s V Duj J J_ oc 

Due to the properties of the functions, the main contribution in the integral comes from the region where 
r oc r s and, therefore, in Eq. (C6) we have used r + r w r s and we defined there = "^J H . In the limit 
of small (2; — > we have y cx Co 1 / 2 — > and, therefore, we can write the integral in the form: 



— " 
<5w 



— / dye- y y— \ dxe «s , (C7) 



or 

~ Js /.00 /■ A „ J- 



^e 96 ^ / dye-y / dx^) 2 ^ ^ e 96 ^F. (C8) 

OOJ Jo J 00J 



The equation for the intercept of the resulting pole is : 



UJ — Us 

or 



As J & = 0, (C9) 



oj — ujs \ 96wr| 

This leads to the solution: 



i-tA-(i + ^)=o. (cio) 
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96r|(£ s + A s ) ' 



u = ^s-h = A S + u s + n<8 » (Cll) 



In another region of ui we assume that the position of the resulting pole u> = &s-h may be close to Cjs , 
As-h — &s + Sw j where for Slo there are two asymptotic limits. The first is the small Slo , ^ << 1 , 
and the second for the 5u> of the order of lo,s , 4£- ~ 1 . We solve Eq. (49)by using the Green function 
given by Eq. (C3) . Due to the form of functions <p(r) the main contribution in integral of Eq. (45) comes 
from these r' w r" w r s . Therefore, we perform the integration in Eq. (45) only over the <j>{r) functions, 
neglecting the term proportional to Sr in the integrand. Using the properties of <f> function, we obtain 

J dr' j dr" <j){r') <j>{r") « 4. (C12) 



This follows from the fact that 



/ 



dr 4>{rf = A / %-e-tl& = 1. 
Is 



When Slo may be small or of the order u>s , ( u> = ujg + Sw ) , i-e. ^ << 1 or ^ ~ 1 , we obtain the 
following equation for resulting pole: 

The solution of this equation for the to = Cos + with Slo <~ tij , up to the order (^) 2 , is : 

2A|r s 



APPENDIX D 



We again consider our variable £: 



At small values of lo we have: 

' m 1 ■ <™> 



We should calculate the contribution of the secondary poles to the second term of the Green function of 
Eq. (71): 



G DIS (r f rl = 9^,rf,n) 

w V*>' x > l _ Asfdr' j dr" 4>(r') G"{r',r") <j>[r") 



We rewrite this equation in the following form: 
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In the limit of small u> we can neglect 1 in the denominator of Eq. (D2): 



u yf ' l >~ V4™w L r fl D \Lu L r H Dj Ai(£;{ = 0) 1 ' 

The secondary pole positions are defined by the zeroes of the Ai(£;( = 0) function in the limit u) — > 0. 

For £ given by the Eq. (Dl) we use the method of steepest descent to estimate this function. Calculated 
saddle points are equal to ± i v 7 ?- They lead to contribution: 

Ai(^;C = 0)^ « ^ sin (^l^| 3/2 )' (° 4 ) 
which together with Eq. (D3) gives: 

C DIS (r r ^ ?j (A(g; Qjr^ [M& Q)r f =r s 

u {rf ' n) ~ (4w L r H D 2 u; 2 y/ 3 sin(§|£|3/2) ■ W 
We obtain for the secondary pole positions the following equation ( see also [13], [25]) : 

sm ( o 77 —7- ) = . ( D6 ) 



.Ag 1 ( 4W V (i4i(£ C))r,=rs C))r,=r 



3 V D 

which leads to 



2 / 4 y uj l r H 

^ = 3\D (D?) 



Expanding sin in Eq. (D5) in the vicinity of ui = u n for each n, we obtain: 

r Dis, x (i»(e;C)) CT (i»(£;C))^ 



(D8) 



Integrating over oj and closing the contour in uj over the positions of the secondary poles we obtain the 
answer for the contribution of a separate secondary pole: 

The total contribution of the secondary poles is the sum over all " n" : 



r<Dis (Y . r ,._ v ((Mt;Q) r ^ rs (Mt-:0)r f =rs)^ |(A)il^ ln _^_ 



(D10) 



28 



[1] E.A. Kuraev, L.N. Lipatov and V.S. Fadin, Sov. Phys. JETP 45 (1977) 199; 

Ia.Ia. Balitsky and L.N. Lipatov, Sov. J. Nucl. Phys. 28 (1978) 822; 

L.N. Lipatov, Sov. Phys. JETP 63 (1986) 904. 
[2] A. Donnachie and P.V. Landshoff, Nucl. Phys. B244 (1984) 322; Nucl. Phys. B267 (1986) 690; Phys. Lett. 

B296 (1992) 227; Z. Phys. C61 (1994) 139. 
[3] D. Kharzeev and E. Levin, Nucl. Phys. B578 (2000) 351, hep-ph/9912216. 
[4] D. Kharzeev, Yu. Kovchegov and E. Levin, Nucl.Phys. A690 (2001) 621, hep-ph/0007182. 

[5] E. V. Shuryak, hep-ph/0101269; Maciek A. Nowak, E. V. Shuryak, and I. Zahed, Phys. Rev. D64 (2001) 034008; 

hep-ph/0012232 and references therein. 
[6] A. B. Kaidalov and Yu. A. Simonov, Phys. Lett. B477 (2000) 163; 

[7] R.A. Janik and R. Peschanski, Nucl.Phys. B586 (2000) 163; B565 (2000) 193, hep-th/9907177; 

R.A. Janik, Phys. Lett. B500 (2001) 118, hep-th/0010069; 

K. Tuchin, Phys. Lett. B497 (2001), hep-ph/0007311. 
[8] M. Genovese, N.N. Nikolaev and B.G. Zakharov, JETP 81 (1995) 625; 

B. Z. Kopeliovich, I. K. Potashnikova, B. Povh and E. Predazzi, Phys. Rev. Lett. 85 (2000) 507. 
[9] B.C. Brower, S.D. Mathur and Chung-I Tan, Nucl.Phys. B587 (2000) 249, hep-ph/0003153. 

[10] L. V. Gribov, E. M. Levin, and M. G. Ryskin, Phys. Rep. 100, 1 (1983). 
[11] A.H. Mueller, J.-W. Qiu, Nucl. Phys. B268, 427 (1986); 

L. McLerran and R. Venugopalan, Phys. Rev. D49, 2233 (1994); D49, 3352 (1994); D50, 2225 (1994). 
[12] J.C. Collins and J. Kwiecinski, Nucl.Phys. B316 (1989) 307. 

[13] E. Levin, Nucl.Phys. B453 (1995) 303; " ORSAY lectures on low x deep inelastic scattering" LPTHE, Orsay 
91/02. 

[14] Yu. V. Kovchegov and A. H. Mueller, Phys. Lett. B439 (1998) 428. 
[15] N. Armesto, J. Bartels and M.A. Braun, Phys. Lett. B442 (1998) 459. 
[16] E. Levin, Nucl.Phys. B445 (1999) 481, hep-ph/9806228. 

[17] M. Ciafaloni, D. Colferai and G.P. Salam, JHEP 9910 (1999) 017 and 0007 (2000) 054. 
[18] M. Ciafaloni, M. Taiuti and A. H. Mueller, Nucl. Phys. B616 (2001) 349, hep-ph/0107009. 
[19] " Regge Theory of Low - pr Hadronic Interaction ", editor L. Caneschi, North-Holland, 1989. 
[20] V. S. Fadin and L. N. Lipatov, Phys. Lett. B429 (1998) 127 and references therein. 
[21] G. Camici and M. Ciafaloni, Phys. Lett. B412 (1997) 396, B430 (1998) 349. 
[22] D.A.Ross, Phys. Lett. B431 (1998) 161. 
[23] G.P.Salam, JHEP 9807 (1998) 019. 

[24] J. Bartels, A. De Roeck, C. Ewerz and H. Lotter, " The 7*7* total cross section and the BFKL Pomeron at 500 
GeV e + e~ - linear collider", hep-ph/971050; 

S.J. Brodsky, F. Hautmann and D. E. Soper, Phys. Rev. D56 (1997) 6957. 
[25] L.P.A.Haakman, O.V.Kancheli and J.K.Koch, Phys. Lett. B391 (1997) 157; Nucl. Phys. B518 (1998) 275. 
[26] C.Lopez and F.J.Yndurain, Phys. Rev. Lett. 44 (1980) 1118; 

C. Lopez and F.J.Yndurain, Nucl.Phys. B171 (1980) 231; 

A. Donnachie and P.V. Landshoff, Phys. Lett. B437 (1981) 408; 

R.K.Ellis, Z.Kcenszt and E.Levin Nucl. Phys. B420 (1991) 517. 
[27] M. Ciafaloni, D. Colferai, G. P. Salam and A. M. Stasto, Phys. Lett. B541 (2002) 314, hep-ph/0204287. 
[28] E. Levin and C. I. Tan, "Heterotic pomeron: A Unified treatment of high-energy hadronic collisions in QCD", 

hep-ph/9302308. 

[29] A. K. Likhoded, N. V. Mokhov and O. P. Yushchenko, "Bare Pomeron In Inclusive Processes", FERMILAB- 
FN-0504, Jan. 1989; P. V. Chliapnikov, A. K. Likhoded and V. A. Uvarov, Phys. Lett. B215 (1988) 417; 
A. B. Kaidalov and K. A. Ter-Martirosian, Sov. J. Nucl. Phys. 39 (1984) 979 A. B. Kaidalov, Phys. Lett. B116 
(1982) 459. 

[30] D. E. Kharzeev, Y. V. Kovchegov and E. Levin, Nucl. Phys. A699 (2002) 745, hcp-ph/0106248. 



29 



